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Outline 



• Difficulty of sizing flaws with thermography 

• Variational Principle in Image Processing 

• Approximate solution for flat bottom holes 

• Application of variational principle to hole 
sizing with multiple time images 

• Discussion of limitations 

• Summary 


Composite "Flat" Bottom Hole 

Specimen 
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Approximate Depths of Flat Bottom Holes 0.05, 0.1, 0.15 cm 
Approximate Diameters 1.27, 0.63, 0.32 cm 



Thermal Response at 0.25 

seconds 
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Profile across center of flat bottom holes 



Thermographic Sizing of Flaws 

• Size of small flaws complicated by in-plane 
thermal diffusion 

• Two common approaches have been taken 

— Point by point time analysis - flaw sizing from earliest 
flaw response (limited in-plane diffusion) 

— Image analysis- deconvolution or other image 
processing technique to remove diffusion effects 

• Alternate approach is simultaneously incorporate 
time and spatial information into analysis 

— Simultaneous application of variational method for 
multiple time responses 





Variational Principle in Image 

Processing 


• Applications: 

- De-noising 

- Segmentation 

- Deblurring 

- Estimation of apparent motion 

- Registration 


• Minimize energy => E(u ) = f F(u,I)dx dy + y G(u) 

- u- desired image 

- I - acquired image 

- E(u) - Energy 

- E(u ; \) - relationship between desired image and acquired image 

- G(u)~ constraint 

- y- constant 
— El - domain 



Variational Principle in Image 

Segmentation 




Chan, T. F. and Vese, L. A., “Active contours without edges," IEEE Transactions on 
Image Processing 10, 266-277 (2001) 



Variational Principle for Sizing 

cr 

Thermographic Flaws Using Single Image 

Minimize energy => E(u ) = 

f n (D(u(x,y),x,y, t) — T(x,y, t)) 2 dx dy + y f n \Vu\dxdy 

- D(x, y, t) -flaw response for flaw with shape u 

- u(x,y) - binary representation of hole ( 0 not hole , 1 hole) 

- T (x, y, t) - measured response 

- E(u) - Energy 

- f \Vu\dx dy- circumference of flaw 

Need fast method for estimating 

- D(u(x, y),x, y, t) 


Approximate Solution for Flat 
Bottom Hole of the Form 




D (u, x, y, t) = 2F ^^_ /( } ^ // 5 u(x,y)e (( * *' )2+ (y y'^^dx'dy' 

D(u,x,y,t) - Thermal response of flat bottom hole 

d - depth of top of flat bottom hole 

S - surface of top of flat bottom hole 

k- thermal diffusivity of material 

K-thermal conductivity 

F-flux at surface 

b-constant 

Note: Blurring of the delamination is function of the time of measurement, not the 
depth of the top of the flat bottom hole 



Thermal Response Slot 0.1 cm 
Below Surface at 0.5 Seconds 



ACS - Approximate Convolution Solution with b=2.1*in-plane diffusivity 
FEM - Finite Element Simulations 
Both size and shape in good agreement 



Amplitude of Thermal Response for 
Different Hole Sizes and Depths 



Time dependence of amplitude in good agreement with approximate solution 



Single Time Slice Analysis 


Euler-Lagrange Equation of single time response 


<5(0O,y))(G * (G * tf(0(x,y))-T(x,y,t) + y div(l^£^)=^ x,y)) 


'\V<f)(x,y)\ 


cLt 


- G- Gaussian convolution 

- H - Approximate Hat Function 

- 5 - Approximate delta function 

- T(x,y,t) - measured flaw response 

. . , d(4>(x,y )) 

- Solving for =0 

dt 


Details in paper 


Relationship of 0(x,y) to 
Thermal Response 


< P(x,y ) 



Thermal response of flat bottom hole 
G * H((p(x,y )) 



Shape of flat bottom hole 


tf(0O,y)) 


Measured Thermal Response 



Difference between G * //(0(x,y)) and Measured Thermal Response 
Determines 

dt 



Results for 0.635cm Diameter 
Flat Bottom Hole at 1 sec 








Measured 

Response 


Estimated 

Response 


Difference 


Estimated 
Image of Top of 
Flat Bottom 
Hole (estimated 
diameter 0.64) 



From Single 4>(x,y ) to Thermal 
Response at Multiple Times 


4>{x, y) 


Thermal 
response 
at 0.5 sec 



Shape of flat bottom hole 
H(4>(x,y )) 


Thermal 
response 
at 0.7 sec 



Thermal 
response 
at 0.9 sec 

A> 


Time Series Analysis 

Euler-Lagrange Equation of a series of time 
response 


S(4>(x,y))(Z n (G * (G * tf(0(x,y))-T(x,y,t n )) + y div( 


V(j){x,y ) 
\V(p(x,y)\ 


)= 


d{(j){x,y )) 
dt 


- T(x,y,t n ) - measured flaw response at t n 
_ solving for 

° dt 


Details in paper 


Amplitude 


Time Response of Amplitude of 
Thermal Response 



Results on All Holes 




Depth(cm) 

Time Series 
Image Analysis 
Depth Estimate 
(cm) 

Diameter(cm) 

Average Single 
Image Analysis 
Diameter 
(cm) 

Time Series 
Image Analysis 
Diameter 
(cm) 
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0.56 
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Discussion 



Most significant improvement is on smallest holes 

Estimation of flaw size depth incorporates flaw size 
information and give good agreement with known 
depths for all flaws 

- Earliest significant responses give best estimates for depth 

Constraint minimizing circumference of flaw tends to 
result in circular flaws 

- Good for this case, however, not so good if flaws are not 
circular 

More work need for delaminations 

- Single delaminations 

- Overlapping delaminations 


Summary 




• Good approximate solution for thermal 
response of flat bottom holes is the 
convolution of the shape of the top of the flat 
bottom hole with a Gaussian 

• Variational method for determining the size 
and depth of flat bottom holes gives a more 
accurate value for the size and depth 


